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ABSTRACT 

The value of accurately knowing the absolute calibration of the polarizing elements in the Near 
Infrared Camera and Multi-Object Spectrometer (NICMOS) becomes especially important when con- 
ducting studies which require measuring degrees of polarization of close to 1% in the near infrared. 
We present a comprehensive study of all previously observed polarimetric standards using the NIC2 
camera on NICMOS. Considering both pre- and post-NICMOS Cooling System observations we find 
variations in the polarimetry consistent with the effects of sub-pixel mis-alignments and the point 
spread function. We also measure non-zero results from unpolarized standards indicating an instru- 
mental polarization of p w 1.2%, 8 w 88°. The lack of polarized and unpolarized standard stars with 
which to perform a comprehensive calibration study means we cannot be confident that the current 
calibration will be effective for a number of recent large NICMOS GO programs. Further observations 
of polarimetric standards are needed in order to fully characterize the behavior of NICMOS at around 
p = 1%. 

Subject headings: instrumentation: polarimeters, methods: data analysis 



1. INTRODUCTION 

While photometry and spectroscopy provide informa- 
tion on spatial distributions, chemical compositions, and 
dynamics, polarimetry adds a valuable extra data dimen- 
sion that can be used to determine the geometry of as- 
trono mical objects and the properties of interstellar par- 
ticles ilAntonucci fc Millerl 119851 iHennev. Raea fe Axon! 
Il99l E asen et all 12003) . Polarized light may intrinsi- 
cally originate from the emission process or it may be 
induced via interactions with a d iffuse medium ijCollettl 
119921 lLandi degl'Innocentil l2002j) . In either case, the 
preferred orientation of the electromagnetic vector, as 
well a quantitative measure of that preference, can be 
determined by placing polarizing elements in the light 
paths of detectors an d deriving the Stokes parameters 
j|Chandrasekhalll96?it) . 

From the ground some care is needed in determin- 
ing the Stokes parameters as the atmosphere intro- 
duces photometric variations with amplitudes similar to 
those expected from the polarized light. Space based 
polarimeters are not exposed to these anomalies and 
can use multiple polarizing filters to gather the re- 
quir ed information. However, multicha nnel polarime- 
ters l)Serkowski. Mathewson fc Fordll 975). which are yet 
to be employed on space based polarimeters, do reduce 
some of the complexities of the data reduction through 
synchronous recovery of the Stokes parameters. In fact, 
with enough care, ground-based polarimeters can mea- 
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sure f ractional polarizations of 10 6 or less ijHough et alJ 
2005). 

Dust frequently enshrouds many interesting astronomi- 
cal objects and inhibits the transmission of optical wave- 
lengths. This effect declines as the wavelength of the 
light increases, thus many studies are carried out in the 
infrared. It is therefore of great advantage to be able 
to carry out polarimetric studies with the Near Infrared 
Camera and Multi-Object Spectrometer (NICMOS) on 
board the Hubble Space Telescope (HST). 

As with many instruments, NICMOS executes pro- 
grams that cover a wide range of required accuracies 
and sensitivities. It is therefore essential to probe the 
polarimetric capabilities of NICMOS, especially in low 
polarization targets, and assess any signatures that may 
be introduced by the instrument itself. We present here 
such an assessment based on the prev i ous ef forts of, for 
examp le. IHines. Schmidt fc Schneiderl lEffll and lHmesl 
(2002). By examining the on-line data archive for NIC- 
MOS we can directly compare the results from polarimet- 
ric standards whose degree and orientation of polariza- 
tion have been determined during other investigations. 
The study is limited to include data obtained with the 
NIC2 (AA = 1.9 — 2.1/xm) camera as the number of 
previous NIC1 (AA = 0.8 — 1.3/xm) studies, excluding 
the calibration data, is low. In addition, these NIC1 
programs generally probe polarizations of ~ 40%, i.e., 
#7264, where high accuracies are not necessarily needed. 
However, future works similar to those carried out here 
may also be prudent for NIC1. 

In § |21 we describe the data and its subsequent reduc- 
tion. The methods used in determining the degree and 
orientation of polarization are outlined in §[3] The results 
are presented in §f^before being discussed and concluded 
in § and § |H1 respectively. 

2. OBSERVATIONS AND DATA REDUCTION 

During the HST second servicing mission in 1997, 
the Faint Object Spectrograph was replaced by NIC- 
MOS. This instrument operates from 0.8 — 2.5/im using 
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TABLE 1 
Details of Observations 



Target 


Type 


RA (J2000) 


Dec (J2000) 


Epoch 


Orientation (°) 


PID 


CHA-DC-F7 


B5V (Pol.) 


10 56 12.91 


-76 35 54.2 


1997-09-13 


-36.48 


7692 




10 56 12.03 


-76 35 52.0 


1998-04-20 


-179.00 


7958 






10 56 12.03 


-76 35 52.0 


1998-07-05 


-106.49 


7692 






10 56 12.13 


-76 35 53.2 


2002-09-02 


-50.83 


9644 






10 56 12.15 


-76 35 53.3 


2003-05-20 


-160.83 


9644 


G191B2B 


DAw (Unpol.) 


05 05 30.80 


+52 49 52.9 


1997-12-24 


-74.32 


7904 


HD64299 


A1V (Unpol.) 


07 52 25.61 


-23 17 50.3 


1997-09-01 


12.71 


7692 


HD283812 


A2V (Pol.) 


04 44 25.16 


+25 31 43.1 


1997-09-28 


14.52 


7692 




04 44 25.16 


+25 31 42.4 


1997-12-04 


-47.48 


7692 


HD331891 


A4III (Unpol.) 


20 12 2.29 


+32 47 45.2 


1997-09-01 


-93.57 


7692 




20 12 2.08 


+32 47 42.2 


1998-04-26 


32.55 


7958 






20 12 2.11 


+32 47 43.5 


2002-09-09 


-109.63 


9644 






20 12 2.11 


+32 47 43.5 


2003-06-08 


0.57 


9644 



Note. — Types taken from SIMBAD (http://simbad.u-strasbg.fr/Simbadl, all other data are 
as they appear in the HST archive I http://archive.stsd. edu/h sf^ T (t-'ol.j and (Unpol.) after the 
stellar type refers to polarized and unpolarized standards respectively. 



three independent cameras. The setup allows for broad, 
medium and narrow band imaging, coronographic imag- 
ing, broad-band imaging polarimetry, and grism spec- 
troscopy. NIC2, the camera focused on in this study, is a 
256 x 256 HgCdTe array with a (y/075 pixel scale giving 
a 19"2xl9"3 field of view <|Thompson et alJll99c| . After 
installation on HST the dewar was allowed to be brought 
up to operating temperature. Unfortunately, this tem- 
perature was not tested on the ground and the subse- 
quent ice expansion provided enough dewar deformation 
to allow contact between one of the optical baffles and 
the vapor cooled shield. The resulting heat sink not only 
altered the foci of the three cameras, but also depleted 
the cryogen by the beginning of 1999 (two and a half 
years ahead of schedule). However, a mechanical cry- 
ocooler, the NICMOS Cooling System (NCS), which was 
successfully installed in March 2002, has since brought 
NICMOS back to life and restored the infrared capabili- 
ties of HST. While the ins trument per formance is compa- 
rable pre- and post-NCS l)Hinesll2002|) . it must be noted 
that NICMOS now displays a different set of characteris- 
tics. The reason for these changes are at present unclear, 
however, for individual observations, these effects can be 
incorporated as changes in the transmission efficiencies 
of polarizing elements. One must also be aware that dif- 
ferent quadrants of the array may exhibit different bias 
levels introduced by resetting the array, plus additional 
bias offsets that ar e related to oper ating temperature and 
main bus voltage l)Bergeronll2005j) . While this "pedestal 
effect" has been reduced by modified flight software, it 
has not been entirely eliminated and may lead to incon- 
sistent polarization measurements. More specific details 
can b e obtained from th e NICMOS Instrument Hand- 
book l|Schultz et alJl2005^ . 

A comprehensive search of the NICMOS data archive 
has been performed with the aim of analyzing all ob- 
served polarimetric standard stars. There are also a num- 
ber of extended sources (e.g. CRL 2688 - the Egg Neb- 
ula) available for use in a calibration study. However, the 
ground truth for such objects is useless as the measured 
polarizations are dependent on resolution. They also 
generally exhibit a high degree of polarization and are of- 



ten variable. Consequently, such objects will not be con- 
ducive to determining the polarimetric accuracies for the 
low polarization targets (< 5%) with which we are con- 
cerned. A prerequisite for the standards to be included 
in the sample is independent polarization measures from 
well calibrated and characterized polarimeters. The po- 
larimetric standards found, and associated details, are 
presented in Table ^ In addition to the data in Ta- 
ble ^ there are also pre-NCS data for two unpolarized 
standards, HD10 700 and HD30652 <|Lerov fc Le Borgnd 
1989: ILerovll993|) . Unfortunately, the program for which 
these observation were made (#7614) called for the in- 
vestigation of circumstellar structures. The actual stars 
are therefore hopelessly overexposed for accurate study 
here. 

The data and calibration files were retrieved and 
passed through the latest versions of the calnica and cal- 
nicb pipelines where appropriate. Each non-destructive 
readout of individual observations was examined. The 
total counts in each mult iac cum frame were inspected 
as a function of the total exposure time. Deviations of 
this curve of growth from linear, especially at very early 
times, were deemed to be evidence for the existence of 
signal persistence. In the cases where the telescope was 
offset between exposures, the areas of the array exposed 
to source photons in the previous visit were also exam- 
ined for persistence. The individual exposures were also 
checked for any remaining cosmic rays or hot pixels not 
flagged by the reduction routine. 

The apphot package in IRAF 6 was used in order to per- 
form aperture photometry on each polarized standard. 
Circular apertures with radii from 0.5 to 50.5 pixels (in 
1 pixel intervals) were centroided on each standard. The 
observable point spread function (PSF) was comfortably 
included in the outer aperture. In each frame the sky 
was sampled (and removed from the target photometry) 
using a disk whose inner annulus was defined by the edge 
of the observable PSF, and whose outer annulus was de- 

e IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under co-operative agreement with 
the National Science foundation. 
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fined by the edge of the array. In every case the sky 
aperture was large enough to avoid the irregularities in- 
troduced by the small number statistics associated with 
Poisson noise. 

3. DETERMINING P AND THETA 

The methodologies behind extracting the degree 
and orientation of polarization have been specif- 
ically addressed for NICMOS by many authors 
(Mazzuca, Sparks fc Axonll998: Mazzuca & Hineslll999l 



TABLE 2 

Polarimetric Calibration Parameters 



Spar ks fc Axonlll999_ iHines. Schmidt fc Schneiderll200fl 
Frme^l2002yi . However, we briefly revisit the linear tech- 
nique for the case of three non-ideal polarizers here. 

The reduced, sky subtracted, instrumental counts (in 
DN s _1 space) through each polarizing element are used 
to define an observed intensity vector of the form a = 
[li h]- The degree and orientation of polarization (p 
and 0) are defined by the Stokes parameters which also 
define a vector of the form b = (I,Q, U). The two vec- 
tors a and b are simply related to each other by the linear 
expression [C]b = a, where [C] is a matrix describing the 
characteristics of the k th polarizer, namely the polarizer 
efficiency (j]k), the actual orientation (in radians) of the 
polarizer (<pk), the fraction of the light transmitted in the 
parallel direction (£&), and the fraction of light transmit- 
ted in the perpendicular direction (Ik)- The transmission 
coefficients are then described by; 

1 



X k = ~t k (l + l k ) 



(1) 

F fe =X fe (cos20 fc ) (2) 
Z k = Xkis^k) (3) 
allowing the matrix [C] to be defined as follows: 



C = 



x 2 m Y 2 m z 2 

X 3 r] 3 Y 3 773 Y 3 



(4) 



The linear expression [C]b — a can then be solved for 
V", Q and U using, for example, LU decomposition. We 
can then solve for p and using Equations and 



p = 100% x 



y/Q 2 + U 2 
I 



1 u 

— arctan — 

2 Q 



(5) 



(6) 



In Equation a 360° arctangent function is assumed. 
In addition, the orientation of the frame has to be sub- 
tracted from in order to retrieve the on sky position 
angle. 

We have coded the above process in IDL in order 
to solve for p and given three lists of apertures and 
ph otometric data. This p rocess has also been coded 
bv lMazzuca fc Hinesl l|1999|) into the IDL routine polar- 
ize. pro. 7 Their routine takes the three polarized images 
and produces two dimensional maps of p and 0, whereas 
the code used here produces the measured polarization 
through defined apertures. The two methods are pre- 
cisely consistent when tested on flat (polarized) data, 

7 http:/ /www. stsci.edu/hst/nicmos/tools/polarizc_tools. html 







tfc 


Filter 


**(") 


Vk h -NCS +NCS 


POLOL 

POL120L 

POL240L 


8.84 
131.42 
248.18 


0.7313 0.1552 0.8981 0.8779 
0.6288 0.2279 0.8551 0.8379 
0.8738 0.0673 0.9667 0.9667 


Note. - 


"-" and 


"+" refer to pre- and post-NCS co 



i.e., one value of I for each polarizer across the entire 
array. 

In addition to the aforementioned process, we have also 
included an estimation of the errors associated with p 
and introduced by the yariance in the count ra te. The 
estimation is taken from iSoarks fc Axon! l)1999|) and is 
defined by the variance and covariance in the (Q, U)- 
plane. 

3.1. The Polarimetric Calibration Parameters 

The first estimates of the polarimetric calibration 
param eters were presented by IHines. Schmidt fc Lvtlel 
(1997). Preffight thermal vacuum tests of the polariz- 
ing optics found each element to have a unique polar- 
izing efficiency and position angle offset. The values of 
4>ki Vk, h and tk were estimated by attaching a polarizer 
to the Calibrated InfraRed sourCE (CIRCE). This pro- 
vided uniform illumination across the entire array at a 
known position angle and magnitude. Following the on 
orbit observations of the polarized standard CHA-DC- 
F7 and the unpolarized standard HD331891, values for 
tk were adjusted based on exposures from two separate 
roll angles and from two separate epochs l)Hinesl 11998. 
IHines. Schmidt fc Schneider! 12000^ . Using the pre-NCS 
coefficients from Table |3 values of p for CHA-DC-F7 
were found to be within ~ 0.2% of the ground based 
values. Post-NCS the values of tk were again adjusted 
(0.8774, 8381, 09667; Hines 2002) and use d in such 
studies as Ueta. Murakawa & Meixner (2005). The co- 
efficient matrix used in this study has been constructed 
from the pa rameters in Table 121 using the latest revised 
values of ^ (|Hinesll2005|) . 

3.2. Testing the Routine 

The IDL routine was extensively tested in order to 
check its performance against known (simulated) polar- 
izations. The testing also provided examples of how 
the instrument itself may affect the observed polariza- 
tion (e.g., noise and the PSF). We have generated data, 
mapped onto the NICMOS pixel grid, which simulates 
a star with a Gaussian surface brightness profile (both 
with and without the effects of Poisson noise), a point 
source with a TinyTim PSF ijKrist fc Hookll2004lh and 
the effects of sub-pixel shifts between the pointing of the 
three polarizers. Small spatial shifts between polarizing 
elements are not unexpected and may allow flux gathered 
in one frame to fall into an inter-pixel gap in another. 

The resulting polarization profiles (which all use the 
post-NCS transmission coefficients for NIC2), i.e., the 
polarization as measured through progressively larger 
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Fig. 1. — Results from testing the polarimetry routine. (Left column) From a Gaussian surface brightness profile with (dotted line) and 
without (solid line) Poisson noise. (Right column) Using the PSFs as generated by TinyTim. The solid lines are for perfectly aligned frames 
whereas the dotted lines show the effects of sub-pixel mis-alignment using a Gaussian profile. 



apertures, are presented in Figure It can be seen 
that the routine behaves as one would predict. Generally 
there are variations plus larger errors on the small scale, 
and constant results for large apertures. The addition 
of the TinyTim PSF is seen to increase the small scale 
variations and errors. Introducing sub-pixel shifts allows 
the small scale variations to grow. In all noiseless cases 
the profile behavior is stable in apertures larger than 0'.'5. 
The best polarization estimates for real data are there- 
fore likely to come from large apertures enclosing most 
of the flux. 

In all cases, to get appropriate estimates for the effects 
of Poisson noise, we have replicated the signal to noise 
ratios typically seen in the standards (~ 50) as well as 
the expected polarizations (~ 1.2%). 

4. RESULTS 

We now pay close attention to each of the targets listed 
in Table ^ The results and details of previous ground 
based polarimetry observations, as well as the archived 
NICMOS results, are presented individually, summarized 
in Table Eland compared in Figure [3] 

In order to correct for the wavelength depen- 
dence of polarization l|Serkowski. Mathewson & Fordl 
1975) and allow a comparison between these NIC- 
MOS measurements and ground b ased result s , one 
can use the "Serkowski curve" llSerkowskil 11973 
IWilking. Lebofskv fc Riekelllllp) to find p(2.05). This 
notation corresponds to the value of p (%) at the 2.05/im 
wavelength of NIC2. However, in both cases where 



this may be necessary, there are direct measurements at 
p(2.04) bv lWhittet etal] il99l . 

We have considered photometry from data which have 
also had the model PSFs, as generated by TinyTim, re- 
moved using the Lucy-Richardson deconvolution. In the 
test cases (CHA-DC-F7 and HD331891) a change of less 
than 0.1% in p was seen inside the PSF. Using large 
apertures there was no difference as all of the flux was 
enclosed without the need for deconvolution. However, 
the advantage of de-convolving the PSF comes when at- 
tempting to visually identify remaining hot pixels and 
cosmic rays. The extra flux that may be introduced by 
such pixels hiding in the artifacts of the PSF will have a 
detrimental effect on the photometry and therefore the 
measured polarization. The residuals between the PSF 
deconvolved frame and frame containing the PSF were 
therefore carefully examined for further evidence of hot 
pixels and cosmic rays. Two dimensional interpolation is 
used to replace any such features. All of the remaining 
polarimetry was carried out with PSF deconvolution. 

In this study, where we can see the polarimetric behav- 
ior with radius, it is not obvious which aperture will be 
consistent with previous studies; photometric aperture 
sizes have not been published. Here we choose to inspect 
each frame in order to report the measured polarization 
in an aperture that includes all of the observed flux dis- 
tributed across the array by the PSF. The aperture size 
is also reported in Table |3| 
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TABLE 3 

Results of This and Previous Studies 



Target 






Previous 








This Study 




Ppre(%) 


^pre( ) 




Ref. 


Epoch 


R(") 


p(%) 


0(°) 




i i q a. n ni 


126 ± 4 


5.98 




1 QQ7 (1Q 1 ^ 


0.34 




116 ± 7 












1998-04-20 


0.34 


2.22 ± 0.03 


86 ± 4 












1998-07-05 


0.34 


0.95 ±0.03 


116 ±2 












2002-09-02 


0.34 


1.05 ± 0.02 


119 ±5 












2003-05-20 


0.34 


1.19 ±0.03 


111 ±4 


G191B2B 


0.09 ± 0.05 


157 ±1 


n/a 


2,3 


1997-12-24 


0.26 


1.2 ±0.4 


90 ±5 


HD64299 


0.15 ± 0.03 




n/a 


2 


1997-09-01 


0.34 


1.46 ± 0.06 


67 ±4 


HD283812 


1.31 ±0.07 


35 ± 2 


6.29 


1 


1997-09-28 


0.79 


3.11 ±0.08 


121 ±3 












1997-12-04 


0.79 


1.40 ± 0.03 


35 ±4 


HD331891 


0.04 ±0.03 


79 ±1 


n/a 


2,3 


1997-09-01 


0.34 


1.32 ±0.02 


106 ±4 










1998-04-26 


0.34 


1.19 ±0.02 


97 ±5 












2002-09-09 


0.34 


1.41 ±0.02 


90 ±3 












2003-06-08 


0.34 


0.77 ±0.02 


79 ±3 



Note. — References: 1. IWhittet etai] 1 119921) . 2. ITurnshek et alj 119901) . 3. ISchmidt. Elston fc LupTel 1 119921) . The values of p pre and 
9 pre are those taken from the literature after the p(X m ax) correction. R(") refers to the aperture radius through which the observation 
was made. 



4.1. CHA-DC-F7 

The polarimetric p roperties of CHA-DC- F7 were orig- 
inally reported by IWhittet et al.1 l|1992() before be- 
ing us ed in the previous NICMOS p olarimetry stud- 
ies of Mines. Schmidt fe Lvtld (11997^ . IHinesI (119981) . 
Hines. Schmidt fc Schneiderl (|2000l) ~and IHinesI (|2002ji . 
Whittet et alJ lll992|1 report, from ground based mea- 
surements, p{\ m ax) = 5.98% at 0.55/im with an 
instrumental polarization of 0.03%, and p(2.04) = 
1.19% ± 0.01% at 9 = 126° ± 4° which we use in 
this comparison. The NICMOS polarimetry studies 
l)Hines. Schmidt fc Schneiderl |2000|) report p(2.05) = 
0.97%±0.2% at 9 = 119°±6° andp(2.05) = 1.00%±0.2% 
at 9 = 119° ± 6° for the 1997 and 1998 epochs respec- 
tively. 

Exposures of 13.95 seconds were made, at every epoch, 
through each of the three polarizing elements. With the 
exception of the 1997 and July 1998 epochs all observa- 
tions used a four point dither pattern, each frame be- 
ing offset by 2'.' 3. All data in this study are consistent 
with previous polarization estimates apart from the April 
1998 epoch. We note that of all the data for CHA-DC- 
F7 these observations are closest (~ 50 minutes) to a 
South Atlantic Anomaly passage and are the only pre- 
NCS observation to be dithered. While the number of 
cosmic ray hits does not appear to be higher than in 
the other epochs, we have noticed large differences in 
the photometry from different parts of the array. Fig- 
ure |3 demonstrates the difference between photometry 
obtained from the individual non-destructive readouts 
from the dithered and non-dithered data. It can be seen 
that there is a larger spread in the photometry from 
the dithered data, which in turn could lead to a greater 
amount of polarization due to the greater flux differences 
between the polarizers. This effect is not present in post- 
NCS observations and suggests that the pedestal effect 
is non-negligible in pre-NCS dithered polarimetry. It is 
well known that the NCS produces much more stable 
temperatures for th e array ijSchultz. Rove fc Sosevll2"003l: 
lArribas et al.12 005). which is why the pedestal effects are 
much less dramatic than when NICMOS was cooled with 



nitrogen ice. We also note that there are late variations 
from linear (~ 10 seconds onward) in the photometric 
curves of growth. This may be indicative of saturation, 
or at least non-linearity, but we do not find any evidence 
for this in the Data Quality (DQ) data extensions. DQ 
values of 3072 are reported which correspond to "Pixel 
containing source" (1024) and "Pixel has signal in the 
0th read" (2048), not "Saturated pixel" (64) or "Poor or 
uncertain Linearity correction" (2) 8 . Following on from 
this we also do not see any non-linearity early in the 
curves of growth. This indicates that persistence is in- 
significant. 

4.2. G191B2B 

ISchmidt. Elston fc Luniel p)^ report G191B2B to be 
an unpolarized standard with p(0.55) = 0.09% ± 0.05% 
and 9 = 157° from observations using the "Two-Holer" 
polarimeter. The instrumental polarization is estimated 
to be 0.05-0.10%. 

Each NICMOS observation has an exposure time of 
23.97s and has been dithered using a three point pat- 
tern. Compared to the CHA-DC-F7 observations the 
point spacings for these dithers are small. In addition, 
the spread in the photometric curves of growth is less 
than those seen in Figure Nevertheless the results 
from this study do show polarization at a level of 1.2% 
where something much closer to zero is expected. 

4.3. HD64299 

HD64299 is listed as an unpolarized standar d 
0.55) = 0.15% ±0.0 3%) bv lTurnshek"et~aT] lfl99(il) . 
lAgeorges fc Walsbl l)1999j) also use this standard in or- 
der to characterize the instrumental polarization of the 
ADONIS polarimeter (p ins (~ 2.0) w 1.5%). 

No dithering was used during each of the 11.96s expo- 
sures and there is no evidence for persistence or satura- 
tion. However, as with the case of G191B2B, we measure 
a polarization that is inconsistent with zero. 

8 See Chapte r 2 of the NICMOS Data Handbook, 

IMobasher fc Rovd <200l) . 
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Fig. 2. — Example polarimetric curves of growth. Error bars 
have been omitted for clarity but are approximately ±0.03% in p 
and ±4° in 9. 
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Fig. 3. — Comparing non-destructive readouts from dithered and 
un-dithered data for CHA-DC-F7. Inset: An expanded view from 
the end of the exposures highlighting the photometric differences. 



4.4. HD283812 

The polarized properties of HD283812 have been re- 
ported to be p(~ 0.55) = 4% ± 1%,6 = 33.8° 
and p(0.55) = 6.29 % ± 0.05%, 9 = 32° ± 1° by 
ITurnshek et all dT990h andlWhittet et all <|l992) respec- 
tivelv. iWhittet etlfll l)1992j) also measures p(2.04) = 
1.31% ± 0.07% at 9 = 35° ± 2°. As HD283812 is an 
extended source, we use the brightest northern compo- 
nent to begin the centroid for photometry. 

Both epochs of the data for HD283812 had exposure 
times of 5.98s. Dithering was only applied to the Septem- 
ber 1997 epoch data however, where the point spacing 
was large (5'.'0). The dithered photometric curves of 
growth show a spread comparable to the dithered data 
in Figure |21 there is no evidence of persistence. We find 
that the dithered data shows a larger degree and orienta- 
tion of polarization than expected but the non-dithered 
data is entirely consistent. 

4.5. HD331891 

ITurnshek et alJ l)199f)D and iSchmidt. Elston fc Lupiel 
(1992) report HD331891 (BD32+3739) to be an unpo- 
larized standard with p(~ 0.43) = 0.04% ± 0.03%. This 
standard was also used in determining the pre- and post- 
NCS coefficients. 

Whereas all exposures were 5.98s, all but the 1997 
data have been dithered with a 2'.'3 four point spacing. 
There is no evidence of persistence and none of the data 
shows an exaggerated spread in the photometric curves of 
growth. Nevertheless, all epochs show non-zero degrees 
and orientations of polariz ation. These results are con- 
sistent with the findings of lUeta. Murakawa fc Meixneil 
( 2005) who use the same data with the same co-efficients 
to find p(2.05) = 1.4%. They attribute the > 1% polar- 
ization, i.e., greater than the S. 1% instrumenta l polar - 
ization reported bv lHines. Schmidt fc Schneiderl $2000), 
to "systematics in the data reduction procedure." How- 
ever, the estimate of pi„ s < 1% was based upon ground 
based-thermal vacuum tests and not from on-orbit data. 

5. DISCUSSIONS 

As demonstrated by Figure |21 we can see that the simu- 
lated (Figure 0) and real data compare well. In all cases 



the large-scale values of p and 9 are clearly stable. The 
inner variations due to the residual PSF and sub-pixel 
mis-alignment are also seen, but are insignificant outside 
of a radius which approximately corresponds to the first 
airy ring of the PSF (~ 0'.'3), i.e., as long as the pho- 
tometric apertures contain all of the flux spread by the 
PSF, consistent values of p and 9 will be measured. 

It is clear that there are significant differences in po- 
larizations measured at separate epochs for CHA-DC-F7 
and HD283812. It is unlikely that these variations are 
intrinsic to the objects. In both cases it is noted that the 
data are pre-NCS and dithered. It is also unclear why 
such an anomaly is not seen in the pre-NCS dithered data 
for HD331891. However, as suggested in H4.ll this may 
be an artifact from the pedestal effect. There do exist 
several IRAF software packages to remove the pedestal 
effect (pedsky, pedsub), but as the dither allows a check 
of the NICMOS polarization properties in each quadrant, 
we can check to see if this is the reason for the discrepant 
results by performing our routine only on the data from a 
single pointing. Accordingly we have re-reduced the raw 
data files through the calnica pipeline and determined 
the polarization results from each quadrant. The results 
are presented in Table 0] It can be seen that the individ- 
ual pointings give polarization measures consistent with 
un-dithered data, and that the pre-NCS pedestal effect 
can induce variations of Ap w 0.2% and A9 m 15° be- 
tween different quadrants. Post-NCS the pedestal effect 
does not appear to have affected the results, however, 
it is clear that the pedestal effect does produce slightly 
different polarizations in each quadrant of the array. 

We have measured polarization in all of the targets 
which have previously been listed as unpolarized. Such 
findings suggest that there may be a low level of instru- 
mental polarization. From Table |3| it can be seen that, 
on average, this instrumental polarization has a magni- 
tude p w 1.2% at a position angle of 9 « 88°. We now 
demonstrate how such an instrumental polarization may 
be corrected for. 

The fact that polarized standards have been mea- 
sured to be consistent with previous studies (CHA-DC- 
F7, HD283812) allows us to be confident that the re- 
sults for the unpolarized standards are real, and that 
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TABLE 4 
Pre-NCS Pedestal Check 



TABLE 5 

Correcting for Instrumental Polarization 



Quadrant 


CHA-DC-F7 


HD283812 


p(%) 




p(%) 


en 


1 


1.25 ± 0.03 


115 ±6 


1.49 ±0.03 


31 ±2 


2 


1.12 ±0.03 


128 ±7 


1.41 ±0.01 


24 ±4 


3 


1.13 ±0.03 


113 ±7 


1.53 ±0.03 


35 ±4 


1 


1.21 ±0.02 


114 ±6 


1.31 ±0.03 


36 ±2 



Target 



Epoch 



Pcorr(%) 



-C) 



we may be seeing a residual instrumental polarization. 
The low number of observed unpolarized standards and 
the pedestal effect make it difficult to be confident in 
the exact nature of the instrumental polarization, but 
we can nevertheless make an attempt to correct the data 
for this. As the degree and orientation of polarization 
are derived in Stokes (I, Q, U) space this is also where 
we can perform an instrumental correction. It is essen- 
tial that the correction is carried out in the (Q, [7)-planc 
for each standard and not in the celestial coordinate sys- 
tem. The unpolarized standard data from pre- and post- 
NCS observations are considered separately. While the 
(Q, U) quadrants are consistent between the polarized 
standards, the actual values of Q and U are dependent 
on /. The corrections are therefore averaged (weighted 
by the errors in Q and U) across all observed standards 
after being normalized by the intensity, rather than being 
simply carried out in the (Q, C/)-plane. The results from 
the correction are presented in Table [3] It can be seen 
that the dispersion around the expected and previously 
reported results has been reduced. The x 2 for the polar- 
ized standards about the expected values is now ~ 0.01 
as compared to ~ 0.1 without the (Q,U) correction. 

It is also possible to null the polarization in unpolarized 
targets by adjusting the on orbit derived transmission co- 
efficients. As Q = U = in unpolarized sources, tk can be 
derived by considering Ik = IXk and fixing £3 = 0.9667. 
Solving tx and t 2 for p = 0% and 9 = 0° in HD331891 
across all photometric apertures results in the co-efficient 
profiles presented in Figure 0] It can be seen in Figure 0] 
that in the central most aperture the co-efficients derived 
from the 2002 epoch (dashed line) are consistent with 
the co-efficients quoted in Table [S] Similar small scale 
variations to that of the observed polarimetric curves of 
growth are also seen. Taking an average of the derived 
2003 epoch co-efficients, outside of the inner arc-second, 
and carrying through the photometric errors, gives us 
*i = 0.8717 ±0.0005 and t 2 = 0.8341 ± 0.0005. Applying 
these adjusted co-efficients (which are within 1% of the 
original co-efficients) to the 2003 data for CHA-DC-F7 
produces a typical change of Sp ~ 0.1% and 69 m 3° 
over the values of p and 9 derived from the original co- 
efficients. 

Adjusting the transmission co-efficients to null the 
measured polarization in a single unpolarized standard 
may provide a "quick fix" for some data, but it is by 
no means a concrete solution to the underlying problem. 
After all, it is possible there have been no changes in the 
transmission co-efficients. Ideally we would be able to 
characterize these effects by observing several polarized 



CHA-DC-F7 



G191B2B 
HD64299 
HD283812 

HD331891 



1997- 
1998- 
1998- 
2002- 
2003- 
1997- 
1997- 
1997- 
1997- 
1997- 
1998- 
2002- 
2003- 



09-13 

04- 20 
07-05 
09-02 

05- 20 
12-24 
09-01 
09-28 
12-04 
09-01 
04-26 
09-09 

06- 08 



1.24 ±0.03 

1.25 ±0.03* 
1.21 ±0.04 
1.21 ±0.03 
1.12 ±0.04 

0.4 ±0.9 
0.34 ± 0.05 
1.32 ±0.04* 
1.38 ±0.09 
0.36 ±0.05 
0.33 ± 0.04 
0.26 ± 0.04 
0.26 ±0.03 



128 ± 7 
117 ±5 
133 ± 5 
109 ±7 
124 ±6 
-7± 19 
-5 ±17 
33 ± 7 
29 ± 5 
6 ±11 
5 ±13 
-2 ±11 
2 ±15 



Note. — Epochs marked with "*" are the results 
from the 1st quadrant of NIC2 only. 



standards at the cardinal angles of the polarizers (this 
greatly reduces the degrees of freedom in the coefficient 
matrix). This would also allow us to see if there is an 
instrumental polarization that is dependent on the ori- 
entation and level of polarization in the object as well 
as the roll angle of HST. For example, it is possible that 
the sensitivity of the mirrors to polarized light is pro- 
ducing some reflective rctardancc. Unfortunately, there 
is, at present, insufficient calibration data available in 
the archive with which to fully check this possibility. 
There does exist data for two polarized objects, CHA- 
DC-F7 and HD283812, at multiple roll angles (see Ta- 
ble but this statistically insignificant sample inhibits 
us from drawing any conclusions about the nature of an 
observationally dependent instrumental polarization. 

While previous polarimetric studies with NICMOS 
have not needed to be concerned with such character- 
istics due to the high levels of observed polarization 
ljUeta. Murakawa fc Meixnerl2005|) , as the capabilities of 
NICMOS are stretched by more and more ambitious pro- 
grams so will the need to accurately measure low levels 
of polarization. In fact, a number of programs needing 
these levels of accuracy have already been executed. For 
example, #10160 (The nuclear scattering geometry of 
Seyfert galaxies) and #10410 (Anisotropy and obscura- 
tion in the near nuclear regions of powerful radio galax- 
ies) will both require high accuracies as AGN generally 
display p m 1 — 5%. In addition, as the James Webb 
Space Telescope will not be flying with polarimetry op- 
tics, NICMOS will remain the only instrument capable 
of performing such high precision imaging polarimetry 
on faint objects. 

The recent polarization calibratio n for the 
Adva nced Camera for Surveys l|Biretta et all 

2004J iBiretta fc Kozhurina-Plataisl 120041 

iKozhurina-Platais fc Birettal l200i |2005|) presents 
an example for the NICMOS calibration. In order for 
the aforementioned programs to be successful, new data 
on a number of both polarized standards, at various 
degrees and orientations of polarization, and unpolar- 
ized standards should be obtained. Pointings should be 
gridded across the detector in order to investigate field 
dependence. A number of exposures, in each polarizer, 
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Fig. 4. — The values of tk required to null the unpolarizcd 
standard HD331891. The value of to, is held constant at 0.9667. 



at different roll angles should also be employed in order 
to replicate the ground characterization. This will allow 
the derivation of element independent transmission 
co-efficients (from point sources) and will allow tests 
for higher order calibration effects. Such thorough 
observations may also shed some light as to the reason 
for the change in behavior from Cycles 7 and 11. 

6. CONCLUSIONS 

We have conducted polarimctric analyses on all 
(un)polarized standards available in the NICMOS data 
archive for the NIC2 camera. The principle aim was to 
determine the behavior of the instrument at low levels of 
polarization. We have thoroughly tested our routine for 
deriving aperture polarimetry of objects and found it to 
be robust. It has also demonstrated that observed small 
radius variation in the polarimetric curves of growth can 
be attributed to the effects of sub-pixel mis-alignments 
between polarizing elements and the point spread func- 
tion. We have found no evidence for persistence. 



Our findings also indicate that there is a measured po- 
larization of p w 1.2% in unpolarizcd targets which may 
indicate an intrinsic instrumental polarization and the 
need for further tweaking of the on-orbit transmission 
co-efficients. Assuming we are detecting an instrumental 
effect we have corrected data for polarized targets in the 
(Q, £7)-plane and found the dispersion around previous 
(ground based) polarization estimates to have been re- 
duced, but not totally removed. We have also derived 
averaged values of tk that null the unpolarized standard 
HD331891 across all apertures, and found t\ and t 2 to be 
0.8717 and 0.8341, respectively, when fixing £3 = 0.9667. 
In addition, we have attempted to investigate the possi- 
bility of an observationally dependent instrumental po- 
larization, but are inhibited from any conclusions by an 
insignificant sample. 

While such levels of residual intrinsic polarization may 
not hamper studies of highly polarized targets, these lev- 
els will have a detrimental effect on studies attempting to 
measure p < 5%. It is clear that a more comprehensive 
calibration study of NICMOS is critical in order for a 
number of HST programs to be carried out successfully. 
The current post-NCS calibration archive (one polarized 
and one unpolarized standard) is insufficient for this to 
occur effectively. An array of polarized and unpolarized 
standards should be observed, at many roll angles and 
in all quadrants of the instrument, in order for the low 
level polarization characteristics of NICMOS to be prop- 
erly and fully investigated, understood and removed. 
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